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ABSTRACT 


The  electronic  interactions  responsible  for  the  observed  electronically 
induced  phase  transformations  in  several  important  classes  of  materials  are 
being  investigated  theoretically  in  close  collaboration  with  experimental 
programs  at  other  institutions.  These  include  destruction  of  superconductivity 
and  onset  of  a normal  magnetic  state  at  low  temperatures  in  ternary  rare-earth 
superconducting  compounds,  charge  density  waves  in  layered  transition  metal 


dichalcogenides , origin  of  phonon  anomalies  in  Pd  and  Pt  metals,  local  density 


T 


-2 


PROGRESS  ACHIEVED 

The  high  level  of  scientific  progress  already  made  during  this  last  year 
on  the  research  proposed  will  become  apparent  from  the  accomplishments  report- 
ed briefly  in  this  section  and  the  list  of  publications  (attached  as  Appendix) , 
which  have  resulted  from  AFOSR  support. 

A.  Destruction  of  Superconductivity  and  Onset  of  a Normal  Magnetic  State 

\ 

at  Low  Temperatures  in  Ternary  Rare-Earth  Superconducting  Compounds 

Last  year  we  indicated  that  the  search  for  materials  having  higher  super- 
conducting transition  temperatures  and  higher  critical  fields  has  been  exten- 
ded to  the  study  of  ternary  compounds.  Recently,  Matthias  and  collaborators 
at  La  Jolla  had  found  a new  group  of  ternary  compounds  which  show  unusual  tran- 
sitions to  the  superconducting  or  the  ferromagnetic  state.  Their  general  for- 
mula is  MRh^B^,  where  M is  either  a transition  or  rare  earth  (three-or  four- 
valent)  element  and  they  are  superconducting  for  M = Y,  Nd,  Stn,  Th,  Er,  Tm, 
and  Lu  and  ferromagnetic  for  M = Gd,  Tb,  Dy,  and  Ho.  This  is  only  the  second 
true  ternary  system  of  high  T£  superconductors  found  thus  far.  The  first  was 
that  of  the  double  Mo-sulphides  originally  reported  by  Chevrel  et  al.  The  new 
system  is  tetragonal  with  18  atoms/unit  call. 

Matthias  et:  al.  have  reported  most  unusual  behavior  for  the  system  ErRh^B^ 
at  low  temperature.  As  one  cools  down  to  8.9  K,  there  is  a transition  from 
normal  metallic  conductivity  to  superconductivity,  and  as  one  lowers  the  tem- 
perature, there  is  another  transition  at  0.9  K which  returns  the  metal  to  its  nor- 
mal, but  magnetic,  state. 

Following  our  proposal  last  year,  we  have  calculated  the  electronic  band 
structure  of  three  of  these  compounds  ( M ■ Y,  Er,  Ho)  to  determine  a possible 
electronic  origin  of  these  transitions  using  the  self-consistent  LMTO  method 
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of  Andersen  and  the  full  18  atoms /unit  cell.  This  is  perhaps  the  only  method 
which  can  handle  such  a complex  problem  with  a finite  amount  of  computer  time 
and  reasonable  accuracy.  These  calculations  have  been  described  as  represent- 
ing a tour  de  force  on  what  is  the  most  complex  band  structure  problem  ever 
tackled.  They  yield  a good  understanding  of  the  origin  of  the  different  phase 
transitions.  Briefly,  we  find  that  the  superconducting  transition  is  governed 
by  the  Rh  4d  electrons  (through  the  Bardeen-Cooper-Schrieffer  interaction) 
which  show  a peak  in  the  density  of  states  at  the  Fermi  energy  and  the  com- 
peting magnetic/normal  state  transition  is  determined  by  the  rare-earth  5d-4f 
intra-atomic  RKKY  exchange  interaction. 

For  the  alloys  from  Gd  to  Ho,  the  large  effective  spin  moment  leads  to  a 
magnetic  interaction  which  dominates  over  the  superconducting  interaction,  and 
so  the  magnetic  ordering  temperature  Tm  is  greater  than  the  superconducting 

transition  temperature  T#.  In  the  Er  and  Tm  alloys,  with  their  smaller  spin 

moments,  the  "effective"  magnetic  ordering  temperature  the  system  would  have  in 

the  absence  o^  the  c set  of  the  superconducting  state  is  reduced  and  so  I > T 

s in 

The  observed  Tm  in  Er  is  smaller  than  this  effective  ordering  tempera- 
ture because  c c superconducting  state  has  been  achieved,  the  conduction 

electron  suscepribil  is  zero  in  principle — at  least  for  the  Rh  electrons. 
Hence,  for  this  compound,  the  magnetic  interaction  is  largely  ineffective  com- 
pared to  that  in  Ho  (and  earlier  elements  in  the  scries).  However,  at  lower 
temperatures  the  RKKY  magnetic  Interaction  is  still  sufficiently  strong  to  or- 
der magnetically  the  localized  Er  4f  moments  to  produce  a sufficiently  large 
exchange  field  to  then  destroy  the  superconducting  state.  This  results  in  the 

low  observed  temperature  of  the  re-entrant  smgnetically  ordered  state  at 


T-0.9  K 


The  transition  metal  dichalcogenides , MX^*  where  X * S or  Se,  an  im- 


portant class  of  materials  because  of  their  platelet  two-dimensional  structure 
which  causes  them  to  exhibit  strongly  anisotropic  and  anomalous  physical  prop- 
erties, superconductivity,  and  most  recently,  the  existence  of  a transition  to 
the  charge  density  wave  (CDW)  state--a  not  too  well-understood  phenomenon. 

As  Indicated  last  year,  we  have  followed  up  our  successful  studies  of  the 
origin  of  the  CDW  in  TaS2  and  TaSe2  with  a highly  accurate  self-consistent 
study  of  Xi^2  done  with  our  recently  developed  numerical  basis  set  self-consis- 
tent linear  combination  of  atomic  orbitals  discrete  variational  method  (SC-LCAO- 
DVM).  Our  results  showed  unambiguously  that  TiS2  is  an  indirect  gap  (0.25  eV) 
semiconductor  and  thus  resolved  several  controversies  concerning  this  important 
material.  Our  results  showed  good  agreement  with  optical  properties  for  ener- 
gies below  16  eV.  Contrary  to  previous  expectations,  a small  indirect  gap 
(0.2  - 0.3  eV)  was  found  to  occur  at  the  points  L and  M in  the  Brillouin  zone, 
with  a larger  direct  gap  (0.8  eV)  at  r.  To  our  delight,  a recent  comprehensive 
analysis  of  transport,  magnetic  susceptibility.  Hall  effect,  resistivity,  Seebeck, 
etc.,  data  yielded  a gap  of  0.2  - 0.5  eV,  in  agreement  with  our  predictions. 

We  then  called  attention  to  the  need  for  a similarly  accurate  study  of  its 
sister  compound,  TiSe2«  Recent  interest  in  the  nature  and  unusual  prooerties 
of  TiSe_  has  centered  on  the  unusual  2a  by  2c  suoerlattice  which  forms  in  TiSe. 

L O O Z 

at  low  temoerature.  Wilson  has  suggested  that  this  might  be  an  example  of  the 
so-called  excitonic  state  rather  than  a charge  density  wave  state  (CDW)  such  as 
that  found  in  the  Group  V layered  dichalcogenides.  Woo  ejt  a_l.  studied  the  in- 
tensity of  x-rays  diffracted  by  the  super lattice,  magnetic  susceptibility. 
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resistivity,  and  far  infrared  reflectivity,  as  a function  of  temperature  in  non- 
stoichiometric  samples.  Their  results  showed  that  the  transition  observed  at 
* 145  ± 5K  is  second  order,  and  they  concluded  that  it  was  apparently  a nor- 
ma 1-to-commensurate  transition  driven  by  the  Fermi  surface  (FS)  which  has  a high- 
ly two-dimensional  character.  Recent  neutron  studies  on  stoichiometric  samples 
found  Tq  = 202  K involving  transverse  atomic  displacement  with  wave  vector 
q ■ (1/2,  0,  1/2);  this  transition  was  interpreted  as  being  driven  by  electron- 
hole  coupling. 

Since  extrapolation  of  our  results  for  TiS^  to  the  case  of  TiSe^  is  crude 
at  best,  and  since  muffin-tin  or  semiempirical  models  are  not  able  to  resolve 
these  questions,  we  undertook  an  ab,  initio  study  of  the  electronic  band  struc- 
ture of  TiSe2  in  order  to  understand  its  various  measured  properties.  Of  direct 
interest  was  to  compare  with  the  first  angle-resolved  photoemission  studies  of 
Bachrach,  e_t  al. , which  showed  conduction  band  overlap  of  the  valence  band  along 
the  r-M  direction  by  at  least  0.5  eV. 

The  band  structure  was  calculated  in  the  local  density- functional  (LDF)  ap- 
proach, using  our  previously  published  numerical  basis  set  LCAO  discrete  varia- 
tional method.  We  find  that  TiSej  is  a semimetal  with  an  indirect  negative  gap 
(0.20  eV)  between  r and  L;  contrary  to  some  previous  expectations,  no  electron 

pocket  was  found  at.M.  The  dimensions  of  the  electron  pocket  at  L (0.20  LH, 

20 

0.25  LA,  and  0.5  LM)  indicate  the  presence  of  7-8  x 10  conduction  band  car- 

3 

riers/cm  in  the  perfect  crystal.  This  is  in  very  good  agreement  with  the  most 
recent  independently  reported  transport  and  angle-resolved  photoemission  experi- 
ments of  DiSalvo  and  collaborators,  and  Traum  and  collaboratores  at  Bell  Labora- 
tories. 

What  may  be  said  about  the  origin  of  the  CLW?  The  several  mechanisms  proposed 
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involve  one  or  another  form  of  "nesting"  of  electron-electron  or  electron-hole 
Fermi  surfaces.  This  includes  the  "exciton-insultaor"  mechanism  which,  as  em- 
phasized by  Wilson  and  Mahajan,  is  akin  to  the  large  scale  "electron  hole" 
nesting  present  in  Cr  and  responsible  for  its  spin  density  wave.  Because  of 

the  degeneracy  of  the  hole  bands  along  p-A  and  their  near  degeneracy  along  p-M 

and  p-K,  the  volume  of  the  hole  pocket  is  substantially  smaller  than  that  of  the 

singly  degenerate  electron  surface  and  hence  direct  "nesting"  of  these  surfaces 

would  seem  difficult  to  achieve.  It  is  possible,  however,  that  large  peaks 

in  the  generalized  susceptibility,  v{q) , arising  from  characteristic  "volume 

effects"  (i.e.,  two  parallel  electron-hole  bands  crossing  the  Fermi  surface) 

provides  the  overscreening  in  these  strong  electron-phonon  coupled  systems  which 

drives  the  observed  instability.  Plausible  candidates  in  TiSe2  are  the  very  nearly 

parallel  K^-r^  and  the  -H^  bands  which  are  separated  by  a zone  boundary  wave 

vector.  Parallel  bands,  not  necessarily  crossing  E , are  not  peculiar  to  TiSe0 

t z 

alone  but  prevail  in  many  similar  compounds  such  as  Ti$2,  TaS^,  TaSe^  and  WSe^. 
Indeed,  a major  peak  in  y^q)  was  found  in  lT-TaS2  and  lT-TaSe2  to  be  caused 
almost  entirely  by  "volume"  effects.  It  therefore  seems  highly  likely  to  us 
that  this  same  mechanism  drives  the  CDW  instability  in  TiSe2-  Expressed  in  terms 
of  the  bonding  in  the  system,  such  an  instability  can  be  described  as  a partially 
screened  interaction  between  the  electrons  at  L (with  their  wave-function  ampli- 
tudes predominantly  on  the  Ti  site)  and  the  holes  at  p (with  wave— function  ampli- 
tudes predominantly  on  the  Se  site)  with  a zone  boundary  phonon  supplying  the 
coupling  momentum  and  carrying  out  the  structural  change  (e.g.,  shortening  of  the 
Ti-Se  bond).  Finally,  whereas  the  phenomenological  correlation  between  the  ob- 
served instability  temperature  and  the  c/a  ratio  (or  lattice  ionicity)  works 
remarkably  well  for  many  layered  dichalcogenides,  including  TiSe2»  it  predicts 
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an  instability  temperature  in  excess  fo  100  K for  TiS0.  Our  analysis  also  suggests 
an  explanation  of  the  failure  to  observe  such  an  instability  in  TiS^;  although 
both  TiS^  and  TiSe2  fulfill  the  volume  effect  condition,  slight  differences  in 
the  ionicities  between  them,  which  are  not  reflected  in  a dramatic  change  in 
the  c/a  ratio,  are  sufficient  to  place  this  "nesting"  in  a gap  region  in  TiS^ 
and  in  a semimetal lie  region  in  TiSe~. 


C . Origin  of  Phonon  Anomalies  in  Pd  and  Pt  Metals 

The  possible  role  of  electronic  structure  on  observed  phonon  anomalies  in  met- 


als remains  an  important  fundamental  problem  in  understanding  lattice  contributions 
to  phase  transitions. 

In  their  classic  paper,  Miiller  and  Brockhouse  observed  an  anomaly  in  the 
slope  of  the  (0§§)  branch  of  the  dispersion  curves  for  palladium.  The  anomaly 
not  only  extended  over  a broad  range  of  wave  vector  space  but  also  decreased  rapid- 
ly in  amplitude  with  increasing  temperature.  Later,  Dutton  et  al.  also  found  a 
rather  similar  but  not  identical  anomaly  in  the  T^  branch  of  platinum.  On  the 
basis  of  the  known  Fermi  surfaces,  these  anomalies  were  suggested  as  arising  from 
the  generalized  Kohn  effect  due  to  nesting. 

We  have  studied  the  possible  occurrence  of  generalized  Kohn  anomalies  by 
direct  calculation  of  the  generalized  susceptibility,  xOl)*  We  find  that  a peak 
in  xCq)  can  be  directly  related  to  the  observed  anomaly  in  the  (05§)  T^  branch  of 
the  dispersion  curves  of  Pt  and  Pd.  We  calculated  the  band  structure  for  Pd  using 
the  APW  method  and  the  relativistic  AlVf  scheme  for  Pt--both  with  a warped  muffin- 
tin  potential.  The  band  structures  were  previously  found  to  yield  Fermi  surface 


radii,  temperature  dependencies  of  the  static  magnetic  susceptibility,  x(T),  re- 
sistivity, and  a spin  lattice  relaxation,  TjT,  in  very  good  agreement  with  experi- 
ment. In  the  x(<l)  calculations,  we  used  2048  tetrahedra  in  1/48' th  Irreducible 
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BZ  and  the  energy  eigenvalues  as  fitted  to  a Fourier  series  representation.  The 
results  of  x(<l)  for  q along  [OlO],  [llO],  and  [ill]  directions  were  calculated 
for  Pd  and  for  Pt  for  bands  4,  5,  and  6 which  cross  the  Fermi  energy.  The  intra- 
band  parts  of  x(q)  at  q = 0 for  both  metals  are  found  to  agree  with  the  density 
of  states  at  the  Fermi  energy  to  within  0.57..  Results  show  that  the  dominant  con- 
tribution to  X.  arise  from  the  dominant  band  5 whose  "jungle-gym”  FS  has 

intra 

strong  nesting  features;  the  main  peak  for  Pd  occurs  at  the  same  q value  (=0.65^/a) 
for  q along  the  [oqo]  , [q,q,o],  and  [q,q,q]  directions.  The  locus  of  this  main 
peak  is  a square  in  the  (0,0,1)  plane.  The  maximum  of  x^ntra  for  q along  the  [llO] 
and  [ill]  directions  are  237.  and  137.,  respectively,  higher  than  the  value  of  x(q) 
at  q = 0.  For  q along  the  [010]  direction,  the  peak  is,  however,  5.47.  lower  than 
the  value  of  xintra  at  q = 0.  Hence,  while  phonon  anomalies  are  predicted  for  the 
[llO]  and  [ill]  directions,  no  anomaly  is  predicted  for  the  [100]  direction.  The 
predicted  q value  for  the  ilO]  anomaly,  q = 0.65  ™/a  is  dose  to  the  experimental 
value  of  ~0.7  Tr/a.  Although  there  may  be  a hint  of  an  anomaly  at  0.56  [ill]  in 
the  measurements,  a more  detailed  investigation  of  this  region  is  called  for. 

For  platinum,  ^ntra  f°r  <1  along  the  [010],  [llO],  and  [ill]  directions  has 
mean  peaks  which  occur  at  q = o.68  tr/a>  0.75  ^/a,  and  0.85  tt /a,  respectively. 

Here  too,  this  main  peak  comes  from  the  nesting  of  the  jungle-gym  Fermi  surface 
which  is  not,  however,  as  flat  as  that  of  palladium.  Anomalies  are  predicted  (al- 
though weaker  than  in  Pd)  along  [llO]  and  [ill]  but  not  at  [l00].  The  [llO]  anoma- 
ly is  close  to  the  measured  q value  (~0.7-0.8  TT/a).  Also  in  agreement  with  exper- 
iment, we  predict  a weaker  [llO]  anomaly  for  Pt  than  for  Pd,  i.e.,  the  ratio  of 

[xintra  <q  max)  ' Xintra  (q  = o)]  toX intra  (q  = o)  18  23%  for  Pd  versus  127.  for 
Pt.  In  both  Pd  and  Pt,  weaker  anomalies  are  predicted  for  the  lill]  direction  than 
for  the  [llO]  direction. 
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We  have  also  calculated  \(q)  of  both  Pd  and  Pt  for  q along  an  of f -symmetry 
direction,  TW.  These  results  also  show  peaks  from  band  5 which  arise  from  nest- 
ing FS  features  but  which  are  lower  than  \(q  = 0).  These  results  are  consistent 
with  Miiller's  observations  for  the  "off -symmetry"  direction  in  the  (001)  plane 

o e 

making  an  angle  0 = 30  and  is  no  longer  visible  at  0 = 45  . 

D . Local  Density  Formalism  for  Electronic  Bands 

We  have  continued  our  development  and  testing  of  the  local  density  func- 
tional formalism  (LDF)  approach  to  the  determination  of  electronic  band  struc- 
tures. This  is  the  approach  we  have  already  used  successfully  for  the  layered 
dichalcogenides  studies  described  above.  We  have  successfully  completed  a paral- 
lel study  of  the  important  semiconductor  system,  CdS . Long  a favorite  II-VI  com- 
pound for  study,  recent  interest  in  the  electronic  structure  of  CdS  derives  from 
its  use  in  solar-cells  (e.g.,  in  CdS-Cu^S  heterojuctions)  and  the  observation  of 
electron-hole  complexes  for  high  excitation  densities  in  this  material. 

There  now  exists  a vast  literature  comprising  experimental  and  theoretical 
investigations  of  the  electronic  structure  and  related  properties  of  the  impor- 
tant II-VI  compounds.  Of  obvious  importance  for  understanding  and  interpreting 
many  of  the  observed  properties  of  these  materials  is  a detailed  knowledge  of 
their  underlying  electronic  band  structures.  In  this  regard,  recent  optical  and 
photoemission  experiments  have  played  an  important  role  in  elucidating  key 
features  of  the  band  structure  over  an  increasingly  large  range  of  energies 
above  the  fundamental  band  gap. 

On  the  theoretical  side,  pioneering  ab  initio  energy  band  studies  were  car- 
ried out  for  many  of  these  materials  using  the  non-muf fin-tin  seK-consistent 
orthogonalized  plane  wave  (SCOPW)  method.  This  technique  has  been  applied  in  re- 
cent years  to  the  study  of  the  electronic  structure  of  a large  number  of  semi- 


conductors  and  semi -metals,  and  has  resulted  in  a considerable  gain  in  insight 
into  the  systematics  of  both  the  III-V  and  II-VI  compounds  a'  well  as  in  a bet- 
ter understanding  of  the  photoemission,  high  pressure,  and  optical  data  for  these 
materials.  At  that  time,  the  only  popular  alternative  method  was  the  non-self- 
consistent  empirical  pseudopotential  approach,  which  avoided  completely  the 
specification  of  the  microscopic  potential  in  the  solid  (i.e.,  the  nature  of  the 
Coulomb  exchange  and  correlation  fields)  by  parametrizing  empircally  the  first 
few  Fourier  components  of  an  unspecified  crystal  potential.  Although  this  pheno- 
menological approach  has  produced  a great  deal  of  information  on  the  variation 
of  the  gaps  in  II-VI  compounds  and  the  origin  of  the  important  band  structure 
critical  points  in  the  Brillouin  zone,  no  meaningful  comparison  of  the  adequacy 
of  the  SCOPW  approach,  as  well  as  that  of  the  underlying  exchange  models,  could 
be  made. 

More  recent  developments  in  semi-empirical  non-local  (non-self-consistent) 
pseudopotential  theory  have  improved  considerably  the  agreement  between  the  cal- 
culated band  structure  and  x-ray  and  UV  photoemission  data  by  introducing  a few 
additional  disposable  parameters  to  describe  the  geometry  of  the  angular  momentum 
dependent  model  pseudopotential  terms.  Although  a successful  parametrization 
scheme,  little  could  be  learned  regarding  the  microscopic  nature  of  the  elemen- 
tary interactions  which  lead  to  this  success.  Further,  due  to  the  independence 
of  the  parametrized  model  potential  on  its  associated  variational  wavefunction, 
questions  such  as  self-consistency  and  the  nature  of  the  crystal  dissociation 
products  (which  are  well-defined  IDF  atoms  in  both  LCAO  and  SCOPW  models)  have 
remained  large  unresolved  and  have  inhibited  a meaningful  comparison  with  first- 
principle  models  based  on  well-understood  (but  less  successful)  crystal  poten- 
tials. Subsequent  studies  on  the  convergence  properties  of  the  OPW  model. 
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independent  theoretical  studies  on  compounds  made  of  first-row  atoms  (diamond, 
boron  nitride)  and  photoemission  studies  on  the  core  levels  of  compounds  con- 
taining the  4d  elements,  have  revealed  serious  difficulties  with  this  represen- 
tation. These  include:  (i)  errors  of  the  order  of  1-3  eV  in  the  calculated 
p-like  valence  and  conduction  band  states  of  compounds  made  of  first-row  atoms 
due  to  convergence  difficulties  in  the  absence  of  a pseudopotential  cancella- 
tion for  these  l = 1 states,  (ii)  spurious  shifts  of  about  12  eV  of  the  core- 
like 4d  levels  in  the  course  of  the  self-consistent  treatment,  and  (iii)  errors 
in  excess  of  6 eV  in  the  level  ordering  of  the  cation  4d  and  anion  s bands  rel- 
ative to  recent  photoemission  data. 

We  have  carried  out  a self-consistent  all-electron  first-principles  elec- 
tronic band  structure  study  of  cubic  CdS  carried  out  within  the  local  density 


functional  (LDF)  formalism,  using  our  previously  developed  linear  combination 
of  atomic  orbitals  (LCAO)  technique  in  a numerical  basis  set  representation. 

In  the  first  set  of  calculations,  we  employed  the  same  lattice  constant, 
exchange  (only)  potential  and  computational  parameters  as  were  used  by  Stukel 
et  al.  in  their  early  self-consistent  orthogonalized  plane  wave  (SCOPW)  inves- 
tigation so  that  a direct  comparison  of  results  can  be  made  and  the  validity  of 
the  SCOPW  approach  for  covelently  bonded  4d  systems  can  be  assessed.  In  the 
second  set  of  calculations,  the  Stukel  et  al.  computational  restrictions  were 
relaxed,  a more  accurate  lattice  parameter  was  employed  and  the  Kohn-Sham  ex- 
change and  the  Singwi  e£  al.  correlation  potential  used  to  obtain  the  local  den- 
sity formalism  solutions  to  the  band  problem,  including  variation  of  the  band 
structure  and  related  properties  with  pressure  (change  of  lattice  constant). 
Comparison  with  optical  and  x-ray  and  UV  photoemission  experiments  for  excita- 
tions of  both  the  s-p  and  metal  d bands  in  the  5-19  eV  region  indicate  very  good 
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agreement.  The  direct  gap  at  T was,  however,  found  to  be  0.5  eV  (257.)  too  small, 
a discrepancy  similar  to  that  previously  found  in  non-empirical  studies  for  other 
heteropolar  insulators  (e.g.,  Ne  and  LiF).  This  is  traced  to  the  neglect  of  the 
different  orbital  relaxation  at  the  r and  band  edges  and  to  the  non-cancel- 
lation of  the  self-interaction  terms  characteristic  of  the  local  density  poten- 
tial. Simple  atomic  total  energy  models  for  these  effects  are  shown  to  bring 

\ 

this  gap  into  good  agreement  with  experiment.  From  these  studies  we  concluded 
that  a first-principles  (parameter-free)  exchange  and  correlation  LDF  model  de- 
scribes very  well  the  main  electronic  structure  features  of  the  system. 

- E.  Metal-Insulator  Phase  Transition  in  NbCL 

Niobium  dioxide  shows  a structural  phase  transformation  (rutile 
to  body  centered  tetragonal)  around  810°C  accompanied  by  rapid  changes 
in  resistivity  and  magnetic  susceptibility.  The  high  temperature  rutile 
phase  is  metallic,  while  the  distorted  phase  is  found  to  be  semiconducting. 

The  transition  is  less  abrupt  than  in  the  isoelectronic  compound  VO^  and 
has  been  characterized  as  second  order  by  neutron  diffraction  studies. 

It  is  possible  to  prepare  good  single  crystals  of  NbC^  in  contrast  to  the 
first  transition  series  compounds,  so  this  system  provides  a good  testing 
ground  for  comparing  theories  of  electron- phonon  interactions  and  electronic- 
ally driven  transitions  with  experiment.  In  addition,  pseudobinary  compounds 
such  as  Ti^Nb^  ^0^  are  found  useful  in  studying  effects  of  d-electron 
occupancy  on  electrical,  optical,  magnetic  and  structural  properties. 

Neutron  and  x-ray  studies  have  provided  accurate  values  for  atomic 
positions  in  both  phases  of  NbO^.  The  low  temperature  distorted  phase  is 
understood  to  be  essentially  equivalent  to  a dimerization  of  the  metal  atom 
chains  found  along  the  rutile  c-axis,  analogous  to  the  famous  Peierls  dis- 
tortion. In  the  case  of  VO^,  x-ray  diffuse  scattering  measurements  by 


Terauchi  and  Cohen  confirm  the  existence  of  soft  phonon  modes  centered  at 
the  point  R = (%,0,^)  of  the  Brillouin  zone  as  predicted  by  a band  structure 
study.  However,  the  mechanism  of  the  phase  transition  in  NbO^  is  less  clear 
and  neutron  scattering  experiments  have  failed  to  detect  low-frequency  soft- 
mode  phonon  branches  particularly  at  the  wave  vector  P *■  thought  to 

be  associated  with  the  static  lattice  distortion. 

We  have  studied  the  electronic  band  structure  and  the  generalized 

\ 

susceptibility,  x(q)>  in  order  to  assess  the  possible  role  of  electronic 

instabilities  in  the  observed  structural  phase  transformation.  Since,  as  is 
well  known,  first  principles  calculations  on  such  systems  of  rather  low 
symmetry  are  difficult,  some  calibration  against  experiment  must  be  made 
prior  to  estimating  the  validity  of  the  theoretically  determined  electronic 
response  functions.  Thus  it  is  especially  important  for  our  work  that 
Raccah  e£  aK  have  made  extensive  optical  reflectance  measurements  on  single 
crystal  NbO^  and  Ti^Nb^  ^0^  and  have  derived  the  associated  absorption 
function  62^)  • These  data  provide  the  basis  for  a detailed  comparison 
between  theoretical  band  structures  and  experiment. 

For  the  band  studies,  we  used  the  recently  developed  linearized- 
augmented-plane-wave  (IAPW)  method  of  Koelling  and  Arbman,  a method  proven 
to  provide  an  efficient  and  accurate  approach  for  complex  structures. 

Using  a rigid  band  scheme,  the  results  of  our  calculations  are  found  to 
describe  accurately  the  optical  absorption  of  the  rutile  phase  of  NbO^ 
stabilized  by  20  at  % Ti.  Thus  encouraged,  we  determined  theoretically 
X(q)  in  the  constant  matrix  elements  approximation  incl  -ding  both  Fermi 
surface  "nesting"  and  so-called  'volume'  effects.  Unlike  the  case  of 
VO^,  no  Fermi  surface  nesting  features  are  observed  and  x(q)  not  found 
to  have  maxima  in  the  region  of  the  point  P in  the  Brillouin  zone. 
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